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We provide a detailed study of gravitational reheating in quintessential inflation generalizing 
previous analyses only available for the standard case when inflation is followed by an era dominated 
by the energy density of radiation. Quintessential inflation assumes a common origin for inflation and 
■JL- , the dark energy of the Universe. In this scenario reheating can occur through gravitational particle 

• production during the inflation-kination transition. We calculate numerically the amount of the 

^ ^ radiation energy density, and determine the temperature T* at which radiation starts dominating 

■ over kination. The value of T* is controlled by the Hubble parameter Ho during inflation and the 

I transition time At, scaling as [ln{l/HoAt)f^'^ for HoAt < 1 and Hl{HoAty for HoAt > 1. 

^ |. The model-dependent parameter c is found to be around 0.5 in two different parametrizations for 



o 



the transition between inflation and kination. 
PACS numbers: 98.80.Cq 



I. INTRODUCTION 



Q ' Over the last few years strong evidence has been found allowing to conclude that the Universe underwent an 
i-S^ . accelerated expansion with negative pressure for at least twice in its history: in its very early infancy during inflation, 
and at the present day. While inflation was introduced as an elegant solution to the horizon and flatness problem and 
to explain the anisotropy of the cosmic microwave background radiation [l|, today's accelerated expansion was an 
^ ] unexpected discovery which strongly favors the existence of a dark energy component contributing a fraction J^de — 
0.7 to the closure density. It is somehow intriguing that both periods of accelerated expansion can be explained by 
, scalar fields whose potential energy dominates the energy density at some time of its evolution. In fact, while for 
' inflation an "inflaton" field is introduced, the dark energy component in the present Universe can be explained in 
. a dynamical way by modifying the standard cosmology with the introduction of a slowly evolving scalar field called 
\ "quintessence" Compared to a cosmological constant, the latter approach has the nice feature of explaining in 
a natural way why radiation and dark energy provide comparable contributions to the energy budget of the present 
Universe, in spite of having very different time evolutions (the coincidence problem), through "tracking solutions" 
^ of the quintessence field. It is then natural to imagine whether it is possible to unify the two scenarios, identifying 
, inflaton and quintessence with the same field As a consequence of this, various models of quintessential 

inflation have been discussed in the literature [6| . 

In the quintessential inflation scenario two main qualitative properties emerge as follows. 

i) The potential V{(j)) needs to account for the large mismatch between the inflationary plateau at the beginning 
of the (j) field evolution (whose natural value Vb is expected to be within a few orders of magnitude of the Plank scale 
mp = 2.4 X 10^** GeV) and the tiny scale of the quintessential tail Vp eventually accounting for the cosmological 
constant today, about 10~^^^ times smaller. As a consequence of this, a crucial requirement of V{4>) is to have a rapid 
fall at the end of inflation, so that at the end of its slow-roll phase the field </) experiences a strong acceleration as 
it "deep-dives" from Vq toward Vp- Thus the Universe undergoes a period of "kination" expansion when its energy 
density is dominated by the kinetic energy of 0, whose potential energy eventually provides the dark energy at later 
time. 

ii) In this scenario the standard reheating mechanism which is usually assumed to create the initial plasma by 
the decay oid> is not at work. However, the mechanism of gravitational particle production is able to reheat the 
Universe 0, |8[ without introducing extra ingredients in the scenario of quintessential inflation. Although it is much 
less efficient than the usual reheating mechanism, there is no difficulty to accommodate the required cosmology of 
radiation domination. Initially at the gravitational reheating temperature, the tiny amount of radiation produced 
by this mechanism is largely sub-dominant compared to the energy contribution from kination. However, since the 
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energy density of kination is red-shifted away by the cosmological expansion much faster than the radiation density, 
radiation eventually dominates at some lower temperature which, in viable models, must be above the McV scale 
in order to preserve the successes of Big Bang Nucleosynthesis. Here we remark that other reheating mechanisms 
like preheating or curvaton reheating [lo| could be realized typically by introducing additional scalar fields in 
quintessential inflation. 

The purpose of this work is to provide a quantitative analysis for the amount of particle production and deter- 
mine the temperature T* assuming the gravitational reheating mechanism. The actual value of T* has important 
phenomenological consequences. During kination the expansion rate of the Universe is larger compared to the usual 
radiation domination case, and this can lead to various non-standard cosmological scenarios. For instance, if T* is at 
the GcV scale the relic abundance of a thermal cold dark matter candidate can be significantly enhanced compared 



to the canonical prediction IjJ, because its decoupling time from the plasma is anticipated, providing the correct 
amount of Dark Matter for values of the cross sections between the DM particle and ordinary matter sizeably higher 
compared to the standard case. This has interesting phenomenological implications not only for the LHC and other 
future collider experiments but also for astrophysical observations like direct or indirect dark matter detection. Indeed, 
the enhanced dark matter annihilation rate which is hinted by the recent PAMELA [H] or ATIC/PPB-BETS [l3l[l4| 



data, can be made compatible with the thermal dark matter production in the context of quintessential inflation 



The presence of a kination era may also have an impact on the properties of the clcctroweak baryogencsis [l5j or the 
thermal leptogenesis induced by the CP-violating decay of a right-handed neutrino . 

Due to its interesting phenomenological implications, we aim to discuss in detail within the quintessential inflation 
scenario. The original calculation of gravitational reheating has been performed in Ref. Q in the case of the inflation- 
radiation transition, and its results have been widely used also in papers discussing quintessential inflation Q . In our 
discussion we will generalize the analysis of Ref. (3] to the case of the inflation-kination transition by performing full 
numerical calculations of the Bogolubov coefficients as is reviewed in Section 2. In Section 3 we introduce two ways of 
parametrizing the transitions between inflation and radiation or kination, in order to discuss how the particular way 
of modeling them may affect our conclusions. The results of our numerical calculation of the radiation energy density 
are presented in Section 4 and the kination-radiation equality temperature is discussed in Section 5, showing its 
dependence on the inflationary scale and the transition time. We give our conclusions in Section 6. 

II. GRAVITATIONAL PARTICLE CREATION: A SETUP 

We briefly review here the mechanism of gravitational particle creation 17, Est in order to set up the relevant 
equations in a suitable form for numerical integration. The main results of this section are Eqs. (jlOp and pip . In a 
spatially flat background the Robertson- Walker metric is written as: 

ds^ = a^{r]){drf -dx'dx'), (1) 

where 0(77) is the scale factor as a function of the conformal time 77. Recall that the conformal time 77 is related to the 
usual time thy d-q = dtj a. Let us consider a massless scalar field $ minimally coupled to gravity. Spatial translation 
symmetry allows to separate spatial and time dependence in the Fourier-transform of the field: (f>k{x,t) = pk{i:)e^^'^ . 
Factoring out the scale factor: Xk = Pfc/a, one finds the equation of motion of the (conformal) time-dependent 
harmonic oscillator: 

Xfc+^fcWXfe = (2) 
a" 

where ujliri) = + V = . (3) 

a 

Here the prime indicates a derivative with respect to conformal time and k is the comoving momentum. Assuming 
the adiabatic condition V <^k^ m the early past and in the late future (so that ^(77 ±oo)=constant) implies that 
the field operator can be expanded in terms of a complete set of positive-frequency solutions. Denoting the solutions 
asymptotically free in the past and in the future by fj and Fj , respectively, we have 

^ - E («^/^ + «]/;) - E {^^pj + ^^i) ■ (4) 
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In the above equation aj and bj are operators that annihilate the in and out vacua, and are connected by the Bogolubov 
transformations: 



k 

h ^ {(^Jkaj + (3*kaU . (5) 



Assuming that at 77 = —00 no particles are present, the initial vacuum |0)i„, which in the Heisenberg picture is the 
state of the system for all the time, is annihilated by the aj operators: aj|0)i„=0. However the physical number 
operator that counts particles in the out-region is Nf^ = 6^6fe, so that the mean number of particles created into 
mode k is (iV^,) (0|6^6fe|0)i„ = |/3jfeP, which implies that particle creation is proportional to the negative- 
frequency content of out states with the boundary condition that the in states contain only negative frequencies. This 
is tantamount to solving Eq.([2) with the boundary conditions: 



Xkiv -00) 



2k 



with k |fc|. The Wronskian condition: XkXt' ^ XkX'k ~ * keeps the correct normalization of states and implies the 
additional relation: 

KP-|/3fcP = l. (7) 
Once the function ui{i]) is given, setting Xk = rke^'^'' the real and imaginary part of Eq. ([2]) can be separated: 

Tkn'^ + 2r'knk = 0. (8) 
where = 0'^,. The corresponding boundary conditions are 

rij] —00) 



/2k 

r'{r] —00) = 

17fc(77^ -00) ==-fc. (9) 

The last equation of ([8]) is separable and readily solved, yielding = ^tt-- Then the problem reduces to solving 
the single second-order differential equation: 

^fe('7)-^)n. = 0, (10) 

whose solutions can be found by numerical integration given the function uikiv) from Eq.Q. Once the values of 
^ki'H 00) and rj.(7y 00) are found, the Bogolubov coefficient can be extracted using a suitable combination of Xk 
and its derivative x'k ■ Making use of the asymptotic conditions ^ and of Eq. ([7]) , one finds 

= 2k = 2k ■ ^^^^ 



In the above equation for 77 —00 the boundary conditions r = l/v2k, r' ~ consistently imply |/3fc| — 0. Note that 
whenever V <^ k^ (i.e. close to the adiabatic regime) a perturbative approach is possible in the calculation of the 
Bogolubov coefficient (3k , which can be approximated by the Fourier transform of the function V 0| : 

■ r+oo 

^'^"2k ^'''''VMdij. (12) 
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Although we will mainly rely on the numerical solutions of Eq. ((10|) for our discussions, we will also give some examples 
where the numerical solution is compared to Eg. p^ . Given the Bogolubov coefficient one obtains directly the 
final energy density produced: 
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(27ra)3a 



d^kk\f3k\^. (13) 



III. GRAVITATIONAL REHEATING FOR INFLATION-RADIATION/KINATION TRANSITION 

The amount of gravitational particle production can be calculated by solving numerically Eq. (|10[) given the potential 
V = —a /a = in oJkiv) of Eq. ([3]) (here R is the Ricci scalar). Its functional forms arc uniquely defined during 

the inflation and radiation/kination periods, while a certain degree of arbitrariness arises from the modelization of 
the transition period. For this reason in our analysis we will adopt two different parametrizations for the transition. 
One is a polynomial parametrization for a^{r]) following the original paper 01, and the other is a parametrization 
of the equation state w = p/p smoothly connecting the inflation {w = —1) and the radiation {w = 1/3) or kination 
(w ~ 1) regimes with an hyperbolic tangent. 

Let us first consider the "standard" situation where an inflationary epoch is followed by a radiation-dominated 
period, which has been first discussed in Ref. 0- Taking the (arbitrary) boundary condition, a{—H^^) = 1, with the 
conformal time normalized as 77 = —Hq^ at the end of inflation, one gets during inflation a{t) = e^"*: 

aiv) = -7^ = (14) 

where we have set x = Hqti. Recalling that the Ricci scalar is given by i? = 3(1 — 3w)H^, one has R = 12Hq during 
inflation. If the universe is dominated by radiation after inflation, the adiabatic condition is automatically verified 
since R = with w = 1/3, and one has a" = which implies that a needs to be a linear function of rj. Note that, 
assuming an abrupt transition at some time 770: i? = l^i/g and R ~ for i] < rjo and 77 > ?7o respectively, one finds an 
ultra-violet divergence in the energy density of gravitationally produced particles because the discontinuous change 
in the metric produces too many particles in the high frequency modes [7[. This is just an artifact of the simplified 
situation of an abrupt transition, neglecting the fact that the transition takes place during a finite interval in the 
conformal time Arj. It is then convenient to introduce the parameter xq = HqAtj, and to parametrize the function 
a'^iv) = f{x) as: 

{l/x^ for x < —1 (infiation) 

flo + aix + a2X^ + a^x^ for — 1 < a; < xo — 1 (transition) (15) 
bo{x + bi)^ for a; > — 1 (radiation), 

where the 6 parameters a,; and bi can be fixed by imposing continuity of f{x), f{x)' and f{x)" in the two points 
X = —1 and X = Xq — 1. Since the corresponding expressions we find are quite involved and they differ from the ones 
given in we give them in the Appendix. 

Let us now apply such a parametrization to the different case when inflation is followed by a period of kination. 
In this case one has p oc and if cx l/a^, so that R cx and the adiabatic condition, i? ^ is verified only 
asymptotically for r/ ^ cx) and a(r/) 00. Recalling H ~ a! ja^ one finds that a^iji) is a linear function of r/ during 
kination. The equivalent of Eq. (jlSp for a transition between inflation and kination is then given by 

{Xjx^ for a; < — 1 (inflation) 

ao + a\x + 020;^ + a^x^ for — 1 < x < xo — 1 (transition) (16) 
&o + b\x for X > xq — \ (kination), 

in which the 6 parameters and bi can be fixed in the same way as in Eq. psp . In this case the corresponding 
expressions are particularly simple, and given by: 

«o = -i + 6, «i=-|7 + 8, a2 = -^ + 3, 03 = -^ ^^^^ 
60 = 3 + 3x0 — 61 = 2 + 3x0. 




FIG. 1: Equation of state w = p/p (upper panels) and function V — —a" /a from Eq. ([3} (lower panels) as a function of 
xo ~ HoArj, for xo — 0.1 (left), xq — 1 (center) and xq = 2 (right). Both cases when the Universe at the end of inflation is 
dominated by radiation (dashed line) and kination (solid line) are shown for comparison. The parametrizations of Eqs. (|15l 
[T6)l are used. 




FIG. 2: The same as in Fig[Tl but using the parametrization of Eq. (|18p. 



In Fig [T] we plot the equation of state w ~ p/p (upper panel) and the function V{x) from Eq. ^ (lower panel) 
for both cases of radiation domination (dashed line) and kination domination (solid line) for xq = 0.1 (left), xq = 1 
(center), xq = 2 (right), using the parametrizations of Eqs. (|15ll6p . Taking into accoimt the approximate expression 
of the Bogolubov coefficient given in Eq. (fT2|) it is possible from this figure to anticipate that an enhancement of the 
process of particle creation is expected whenever <C 1, and especially for kination compared to radiation. On the 
other hand, at large xq the process is generally suppressed compared to the previous case and is almost independent 
on Xq. Notice however that the parametrization (|15ll6p can be considered as a reasonable choice capturing model- 
independent properties of a generic transition only when xq 1. When xq > 1 the actual behavior of the functions 
w{x) and V{x) will depend on the details of the model-dependent transition. For instance, the peculiar behavior of 
the function 'w{x) for the inflation-kination transition in the third panel of Fig. [T]is to be ascribed to an artifact of 
the parametrization used rather than to the physical properties of a realistic system. 

In order to estimate the model dependence of our results, in the following analysis we will also adopt a smoother 
transition applied to the equation of state: 
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where y = log(a) + const, and j/o parametrizes the duration of the transition in e-foldings, while Wf = 1/3, 1 is the 
asymptotic equation of state at late times for radiation and kination, respectively. With the above paramctrization, 
the relation 

dlogH dlogH 3 

— 1 = -71 + 19) 

ay d log a 2 

can be integrated analytically to yield: 

H = H/Ho = e-i(^+^i)[y+'^{'°s-°^^^w)] ^ (20) 

with Hq ~ H{y = —oo). The relation between y and x is given by dx/dy = l/{aH). Identifying the transition period 
with the interval —yo < y < +yo (u>(±?/o) = ±0.96), we have: 

.0 = r - (21) 

J-y,H{y')aiy') 

The quantities w = p/ p (upper panels) and V{x) (lower panels) in this paramctrization are shown in Fig. [2] for 
xo = 0.1, 1, 2, where V{x) is obtained from the general relation: 

V{x) = -- = - ^ {l-iw{y)). (22) 

In Fig. [5] we have chosen the (arbitrary) normalization of the scale factor a oc and the origin of the conformal time 
variable x in such a way that the minimum of V{x) is for x = —1, and V{x ~ —1) ~ —2, in order to have a direct 
comparison with the previous paramctrization discussed in Fig. [TJ That is, our normalization in terms of the variable 
y is set to satisfy V{ymin) = —2 and x{ymin) = —1 where ymin can be calculated by solving dV/dy\y=y^.^ — 0. From 
Eq. (pn|) the asymptotic behavior of the Hubble constant follows in a straightforward way: 

H{y » +yo) ^ e-i(i+-^> = , (23) 

which gives, for the case of kination, 

H{y:$> +ya) ^ ^ with h = 2a{yrmnfe-^y"^- . (24) 

All the quantities discussed above can be expressed explicitly in terms of yg. In the case of a transition from inflation 
to kination {wf = 1), we get: 



'(9yo + 6)* + (2/o + 6)i 



^-^'Tl°sl,(9,, + 6)i-(.o + 6)lj' ^''^ 

from which explicit forms of other variables can be derived. 

Note that in the limit j/o ~^ a power expansion allows to recover the paramctrization of Eq. (|16p . In this case 
Ci{ymin), H{y„iin) — > 1, but whcu making a direct comparison with the paramctrization of Eq. ()16p the correspondence 
between j/o smd xq is somewhat arbitrary since w = ±1 only when y ±00. However, by choosing the boundaries 
of the transition at y = ±yo = icyo in such a way that ymin — > —yo — —xo/'2 (so that in both parametrizations the 
beginning of the transition coincides asymptotically to the minimum of the potential), we find 

61 ^ 2e3^° ~ 2 + 3x0, (26) 

which corresponds to the expression for bi found in the last line of Eq. (|17p (The above procedure corresponds to 
setting c = log(?/o)/4)). With such normalizations Figs. [T] and [2] can be directly compared. As expected they differ 
when Xq > 1 while they coincide when ^ 1. Notice however that the definition of the conformal time depends on 
the arbitrary normalization of the scale factor a. For this reason in the following section wc will discuss our results 
expressing the duration of the transition xq in term of the physical time t. Of course the physical results will not 
depend on these conventions. 
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As we will see later, bi is one of the important parameters determining the temperature of the radiation component 
in quintessential inflation models. In the limit of an abrupt transition (At — > 0), we have 61 « 2 as discussed above. 
In the limit of a slow transition, 61 behaves like: 

61 cx (HoAty, (27) 

where At is the transition time and the exponent c is a model-dependent parameter. We find c = 0.5 and c « 0.48 in 
the case of the polynomial and tank parametrization, respectively. For the latter, we define At by At = dy/H. 



IV. NUMERICAL RESULTS 

In this section we discuss our results for the amount of radiation (|13p with the Bogolubov coefficient |/3fcp given by 
the asymptotic value of Eq. (|11[) at rj 00. The latter is obtained in a straightforward way by solving numerically the 
second-order differential equation pop with the boundary conditions given in ^ . In Fig. [3] we present the quantity 
k'^ X |/3fcP, which enters the energy integral of Eq. ([TS]), as a function of k for xq = 0.1 in the case of the inflation- 
radiation transition (thin lines) and inflation-kination transition (bold lines). In the left-hand panel the transition 
is parametrized according to Eqs. (|15ll6p . while in the right-hand one Eq. (fT8|) is used. In both figures solid lines 
are the result of numerical integration of Eq. (|10p , while dashed lines are obtained by making use of the perturbative 
approximation given in Eq. (|12p . A couple of comments are in order here. 

Fig.[3]shows an infrared divergence for the inflation-radiation transition, but not for the inflation-kination transition. 
Such a property can be studied analytically in the limit of an instantaneous transition for which we can obtain simple 
solutions for particle creation. We assume that the inflation-radiation transition occurs at x = Horj = — 1 (xq = 0). 
Apparently natural choices of the in-modes are: 



which are the solutions of Eq. (2) during inflation. Now using the out-modes during radiation in Eq. (6), it is easy to 
flnd 

1/3^1' = (29) 

which results in the infrared divergence of the energy density (fT3|) for the modes fc ^ 0. This implies that the in- 
modcs for k ~* correspond to unphysical states since otherwise this large infrared contribution to the energy density 
cannot represent a self-consistent solution to the Einstein's equations. One way forward to solve this problem is to 
find a suitable choice of the initial vacuum state such that the two-point function for the field does not diverge in the 
infrared limit [l9| . Once such choices are made to obtain finite results for fc ^ 0, their contribution to the integration 
of the total energy density becomes negligible. In this paper, we take a practical approach to cut off the k < Hq 
modes to calculate the energy density. A possible alternative approach to solve this problem would be to calculate 
the radiation energy density by using the renormalized energy momentum tensor with a proper regularization scheme 
[l3| , which would not display this infrared behavior during the quasi dc Sitter stage [2o| . 

It is interesting to note that the situation is different for the inflation-kination transition. Again assuming an 
instantaneous transition at a; = —1, one gets V = Hq/{3 + 2x)'^. While we can take the in-modes as in Eq. ([28|) during 
inflation, the out-modes during kination can be written as 



out _ / ^ 

2 V^o 



akH^^^Kz) + PkH^^\t,z)\ , (30) 



f 1 2) 

where k = k/ Hq, z = .t + 3/2 and iJg ' are the Hankel functions. Matching the in-modes and out-modes at a; = — 1 
one readily finds for fc — > 0, 

(31) 

which is consistent with our numerical results in Fig. 3. Therefore, we can conclude that the radiation-kination 
transition may be free from infrared divergences in reasonable models for the transition. 
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FIG. 3: X l/3fcp as a function of k normalized to Hq in the case of inflation-radiation transition (thin lines) and inflation- 
kination transition (bold lines), and for xq = 0.1. Solid curves are the result of numerical integration of Eq. (|10p . while dashed 
lines are obtained by making use of the perturbative approximation given in Eq. (|12|l . In the left-hand panel the transition is 
modeled according to Eqs. (|15I16|) . while in the right-hand one the parametrization of Eq. (|18|l is used. 

Another interesting observation is that the epoch of particle production is confined to a short interval close to the 
minimum of the potential V{x), which, by our convention, is for x = —1 or y = Umin- This can be seen by considering 
the adiabaticity parameter uj' juj^ which becomes much smaller than 1 for the adiabatic regime. In the case of a slow 
transition (x'q 3> 1) the adiabaticity parameter is suppressed at all times except for a; ~ — 1, so particle creation is 
a quasi-instantaneous process even in this case, as in the situation of an abrupt transition with <C 1. For this 
reason when the transition time is large the amount of particle creation becomes a very slowly varying function of xq 
or HoAt, as will be shown below. 

The computation of the energy density is shown in Fig. 31 where the quantity a'^p normalized to Hq, is plotted 
as a function of HoAt in both cases of a transition between inflation and radiation (dashed lines) and inflation and 
kination (solid lines). Notice how, as previously anticipated, the comoving energy density a^p shows a plateau at 
large At, while an enhancement is present at smaller At. In the case of kination, due to the fact that the function 
V has a second peak at positive 77 for small At (see the lower part of the left panel of FiglT]), which is not present 
for radiation due to a vanishing Ricci scalar R, the quantity a'^p is larger by about a factor of two compared to the 
case of radiation, and in both cases, we find the behavior of ultra-violet divergence: a*p oc In(iJoAi) for At as 
discussed previously. However, the amount of radiation produced by gravitational reheating for large At is similar 
in the two cases, since in that case the potential V for kination is strongly suppressed for 77 > (see the lower part 
of the right-hand panel in Fig. [T|) and thus very similar to the case of radiation. From Fig. [4] we obtain the typical 
amount of particle creation parameterized by / = a^p/ H^: 

0-03 < / < 0.08 Inflation Radiation 

0-03</<0.19 Inflation ^ Kination, (32) 

restricting ourselves to HoAt > 10^^. 

V. RADIATION-KINATION EQUALITY 

A crucial element for the phenomenology of quintessential inflation is the epoch when the standard cosmology with 
radiation domination starts. The amount of radiation produced by gravitational reheating is very small, so that at 
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FIG. 4: The energy density per scalar degree of freedom produced by gravitational reheating (normalized to Hq) as a function 
of the duration in physical time of the transition (normalized to 1/Ho). Solid lines are for the inflation-kination transition 
and dashed lines are for inflation-radiation transition. In both cases the upper line is the result of a numerical integration of 
Eq. (|10|l and the lower one is obtained by making use of the approximation of Eq. (|12|l . In the left-hand panel the transition 
is modeled according to Eqs. (|15ll6p . while in the right-hand one the parametrization of Eq. (|18p is used. 



the end of inflation the energy of the Universe is dominated by the kinetic energy of the quintessence field. However, 
the energy density of kination px is red-shifted away much faster than that of radiation as 

Pk oc and pn cx a"**, (33) 

and thus the latter can eventually dominate at some temperature much below the initial (reheat) temperature Tj. 
Let us define the kination-radiation equality temperature by 

PKin) = PRin), (34) 

which should follow the constraint: > 1 MeV, in order not to spoil nucleosynthesis. 

Consider thermalization of the gravitationally produced particles occurring at a/ with the corresponding (reheat) 
temperature Tj. At this moment, the radiation energy density can be expressed as 

I giiHj , 

Pr = — = ^aiTj , 35) 

ttj 30 

from which one determines a/T/ given I and Hq. The specific value of T/ can be explicitly calculated if the interaction 
rates of particles are given in a particular model 2l|. For instance, assuming a typical interaction rate of T ~ a^T 



with a certain coupling strength a, and equating it with the Hubble parameter H = 6i/2a^ during kination at a/ or 
T/, we find 

, \ — / 3X9 
/QOXs I OT-i \ 

Ho, (36) 




which can be much smaller than the Gibbons-Hawking temperature i7o/27r [22| . Irrespectively of the specific values 
of a/ or T/, we can determine the temperature when the radiation domination starts. Using Eq. (|35p and the 
relation of iso-entropic expansion, gj a^Tf = a^T^, one obtains 



gi_y /30/ 



Ho- (37) 
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FIG. 5; The combination entering into the kination-radiation equahty temperature T, given in Ea. (l39p . as a function 

of the transition time At/Hg^ expressed in absolute time. SoUd hnes are for the transition given in Eq. (|18|l . and dashed hnes 
are for the transition in Eq. (|16|l . In both cases the upper line is the result of a numerical integration of Eq. (|10p and the lower 
one is obtained by making use of the approximation of Eq. (|12|) . 



Since the energy densities of radiation p|j and kination p*j^ at a* or are 



pj, = ^mlHl-^, (38) 



the conditions ([M)) and ([57]) enable us to get 

T.^^^^M. (39) 

In the above equation bi is given by Eq. (|17p or (|24p . depending on which parametrization for the transition is adopted. 
Moreover, notice that both quantities / and 61 depend on the arbitrary normalization of the scale factor, while, as 
expected, in Eq. ((5^ the combination I^/^ /bi entering the physical quantity or Tj does not. 

In Fig. Othc combination I^/^ /bi is plotted in terms of the transition time /St/H^^ . In the figure the solid lines 
show I^/^ /bi when the transition is modeled according to Eq. p^ . while in the dashed lines the transition is modeled 
according to Eqs. (|15ll6p . In both cases the upper line is the result of a numerical integration of Eq. PT7)1 and the 
lower one is obtained by making use of the approximation of Eq. p^ . Note that the quantity I'^^^/bi shows the 
scaling behavior: 

jl r 0.06 [In f for At ^ 

61 \ °iUHo^t)-° '^'' for At ^00, 

where the two values in the second line correspond to the transitions given in Eqs. (16) and (18), respectively. From 
this, we find the following numerical value of T* : 

T ^J^mh-^^l^^WA-f^ for At^O 
* Vff^mp I 0;03g(ij(jAi)=° '« for At ^ 00, 

which depends two important model parameters: the inflation scale Hq and the transition time At. 
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VI. CONCLUSIONS 



In this paper we have provided a detailed study of gravitational reheating in quintessential inflation generalizing 
previous analyses only available for the standard case when inflation is followed by an era dominated by the energy 
density of radiation. In the quintessential inflation scenario both inflation and dark energy are caused by the same 
scalar field: initially the inflaton/quintessence potential energy drives inflation, which is then followed by a kination 
period when the energy density of the Universe is dominated by the kinetic energy of the field; this kination period 
is eventually ended by radiation domination, while later on the potential energy of the quintessence field prevails 
again, providing the dark energy observed today. In quintessential inflation, among several ways to produce radiation 
and thus reheat the Universe, gravitational particle creation is the minimal scenario being able to provide a sufficient 
amount of radiation in spite of a much lower efficiency compared to other mechanisms, thanks to the fact that the 
kination energy density is red-shifted at a much faster rate compared to radiation. In this scenario, the kination- 
radiation equality temperature is an important quantity which has implications for various cosmological events 
like baryogenesis, dark matter decoupling and nucleosynthesis, etc. 

Motivated by this, we performed a numerical analysis of Bobolubov transformations to calculate the radiation energy 
density produced gravitationally, and determine the dependence of T* on the inflation scale Hq and the transition time 
At of inflation to kination. To study how our conclusions depend on the details in modelling the transition period we 
considered two different parameterizations: a polynomial expansion of 0^(77) and a parametrization of the equation 
of state making use of an tangent hyperbolic function. Our main results are summarized in Eq. (|4ip . In both cases 
we obtain numerically similar results, both in the small and in the large At limit which may hold for generic models 
with a reasonable transition behavior. It is also interesting to observe that the inflation-kination transition docs not 
show an infrared divergence for natural choices of the in-modes which represented unphysical vacuum states for small 
k in the case of the inflation-radiation transition. Nevertheless the problem of such an infrared divergence would be 
removed in a consistent way by using the renormalized energy momentum tensor, which could provide an independent 
study of the gravitational reheating process to confirm our results. We leave this issue to a future investigation. 
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APPENDIX A 



We give here the explicit expressions for the parameters and bi introduced in Eq. (|15p for the case of the 
inflation-radiation transition: 



ai 
02 
as 
bo 
bi 



- 3ao - 10 
: 3ao — 15 
flo - 6 

3[ao(l+t) -5-6t] 

26o ' 



where cq is the physical solution of the quadratic equation: 



(AI) 



Cial + 6*200 + C3 = 0, 



(A2) 



chosen in such a way that f{x) = a (a;) is always positive. In (IA2p the constants Ci are given by 



C2 = K^Ks + K2K7 - K4K5 - K^Ke 
C3 = K2KS - KeK4, 



(A3) 
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where 



Ki 


= 3 + 6t + 


K2 


= -10 - 30* - 18*2 


K3 


= 6 + 6* 


< 


= -30-15* 




= 2 + 6* + 6*2 + 2*3 


Ke 


= -20* - 30*2 _ i2t^ 


K- 


= 3 + 6* + 3*2 


Us 


- 10 - 30* - 18*2, 



(A4) 



and * = xo — 1. The above equations differ from Eqs. (18,19) in Ref. [7[, which do not yield continuity of f'{x) and 
f"{x) in X = a;o — 1, as we checked by numerical inspection. 
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